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Abstract 

The in-medium chiral condensate is studied with a new approach which has the 
advantage of no need for extra assumptions on the current mass derivatives of 
model parameters. It is shown that the pion-nucleon sigma term is 9/2 times 
the average current mass of light quarks, if quark confinement is linear. Consid- 
ering both perturbative and non-perturbative interactions, the chiral condensate 
decreases monotonously with increasing densities, approaching to zero at about 4 
fm~^. 

Key words: chiral condensate, sigma term, in-medium 
PACS: 21.65.+f, 24.85.+p, 12.38.-t, ll.30.Rd 



Email address: gxpeng@lns.niit.edu, gxpeng@ihep.ac.cn (G. X. Peng). 
Preprint submitted to Elsevier Science MIT-CTP 3516 



The calculation of in-mcdium quark condensates is of crucial importance to the 
chiral property of quantum chronmodynamics (QCD) [1]. A popular method 
for this is the Feynman-Helmann theorem which gives [2] 
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where E is the energy density above the vacuum, mo is the current mass 
of u/d quarks, and |gg|o = \{(iq)o\ = — (^Q')o ~ (225 MeV)"^. Eq. (1) can 
be schematically derived as such. From Feynman-Helmann theorem one has 



(^(A)||-if(A)|^(A)) = |-(^(A)|if(A)|^(A)) for the hamiltonian H{\) with 



the eigenstate |*(A)). Write the QCD hamiltonian density for the two flavor 
symmetric case as i^qco = H' + 2mQqq where the second term breaks the 
chiral symmetry explicitly while the first term is the remaining part which has 
nothing to do with the current mass tjiq. Then making the substitution A ttiq 



and H Jd^xHc^cD gives 2(^(mo)|gg|^(mo)) = a^(^(mo)|i/QCD|^(mo)). 



Here the integration over space is canceled due to the uniformity of the system. 
The derivative has been changed to partial derivative because the system 
energy may depend on other independent quantities for example the density 
rib. Appply the equality just obtained to the state = |nb) and |^') = |0), 
and then taking the difference lead to Eq. (1) naturally. 

Because no one can exactly solve QCD presently, the energy density E in Eq. 
(1) is given with some model parameters. The main difficulty in this formula is 
that one has to know the current mass derivatives of model parameters which 
are, except for a few special cases, usually not available. 

Recently, another approach has been proposed [3,4] with the advantage of no 
need for assumptions on the current mass derivatives of model parameters. 
However, how the effective Fermi momentum is connected to density was not 
investigated clearly. In Ref. [3], it is merely boosted by a simple factor of 
2 while not boosted in the section 2 of Ref. [4], compared with the non- 
interacting case. In this paper, it will be generally proved that the relation 
between the effective Fermi momentum and density is determined by general 
principles of thermodynamics. At lower densities, a simple relation among 
the pion-nucleon sigma term, the quark current mass, and quark confinement 
is found, which means that, if the confinement is linear, the pion-nucleon 
sigma term should be 9/2 times the average current mass of light quarks. 
Full density behavior of the in-medium chiral condensate is calculated by 
considering both perturbative and non-perturbative interaction effects, which 
gives a zero condensate at about 4 fm~^. 

Let's outline the key points of the new approach in Refs. [3,4] for the two 
flavor case. Ignoring terms breaking flavor symmetry, the corresponding QCD 
hamiltonian density can be schematically written as i^qco = -f^k+2mo^?+i?i. 
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where i^k is the kinetic term, Hi is the interacting part, and niQ is the average 
current mass of u and d quarks. Define a new operator i/cqv = -f^k + 2mqq 
with m being an equivalent mass to be determined by the requirement that it 

meets (ifeqv)nb - {Heqv)o = (-f^QCD)nb " (-f^QCD)o- 

Since particles considered are uniformly distributed, or in other words, rib 
has nothing to do with space coordinates, one can write (\l'|m(nb)^g'|^') = 
'm{n\y){'^\qq\'if). This equality is especially obvious if it is considered in terms 
of quantum mechanics: |\E') is a wave function with arguments rih and coor- 
dinates, the expectation value is nothing but an integration with respect to 
the coordinates. Therefore, if nb does not depend on coordinates, the func- 
tion m(nb) is also a coordinate-independent c-number, and can naturally 
be taken out of the integration. However, if rzb is local, the case becomes 
much more complicated and it will not be considered here. Now substituting 
the expressions of Hqqd and ifgqv into the equation (ifeqv)nb ~ (-f^eqv)o — 

(-f^QCD)nb - (-f^QCD)o ^ads tO 

m = mo + 7T77^\ T^TT = "^o + m- (2) 

Therefore, considering quarks as a free system, i.e., without interactions, while 
keeping the system energy unchanged, quarks should acquire an equivalent 
mass of the form shown in Eq. (2). Prom this equation one can see that the 
equivalent mass m includes two parts: one is the original mass or current mass 
mo, the other one is the interacting part mj. Obviously the equivalent mass is a 
function of both the quark current mass and the density. At finite temperature, 
it depends also on the temperature as well. Because the hamiltonian density 
Heqv has the same form as that of a system of free particles with a density 
dependent mass m, the corresponding dispersion relation e = y/p"^ + m^ is also 
density dependent. The energy density of quark system can then be expressed 
as 



Pf 3 







where (? = 12 is the degeneracy factor, and the function F is defined to be 
F{x) = ^[x^/x^TTi2x'^ + 1) - sh~'^{x)]/x^ with sh"^(,T) = ln(a; + ^/x^Tl). 
The effective Fermi momentum pf is related to the chemical potential /j, by 



Pi = J — m? 01 ^ = Jpj + w?. (4) 



The physical meaning of the equivalent mass is that quarks should have the 
equivalent mass if the system is free but with unchanged energy. It is not 
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difficult to prove that such an equivalent mass always exists. In principle, if 
one obtains the energy density E from some realistic models or even from 
QCD in the future, the equivalent mass m can be obtained by solving Eq. 
(3). Prom the point of view of thermodynamics, one can alway choose to ex- 
press some quantity as an ideal gas form with density and/or temperature 
dependent particle mass(es) while all other quantities determined by thermo- 
dynamic principles. For example, in Refs. [5,6], the pressure of gluon plasma is 
expressed as the same form with that of an ideal gas while the gluon 'thermo 
mass' m(T) is determined by fitting to lattice data. In Ref. [7], the entropy is 
expressed as the same form as that of an ideal gas. Here the energy has been 
expressed as the same form as that of an ideal quark gas. However, how the 
Fermi momentum pf is connected to density should be studied carefully. If no 
interactions are accounted or a specific duality principle is assumed, the Fermi 
momentum will be connected to density simply by pf = [{lS/g)7i'^nb]^^^ . Here 
quark interactions are included within the equivalent mass. So its relation to 
density should be determined by general thermodynamic principles. It can be 
shown that the baryon number density is related to the Fermi momentum 
Pi by 



_ 9P^ 9 



m sh (pf/m) 

Pi - 



mdm. (5) 



To prove this, let's write the fundamental thermodynamic relation d{VE) — 
Td{VS) — PdV + 3fid{Vn\,) which is the combination of the first and second 
laws of thermodynamics. Here P is the pressure, S is the entropy density, and 
V is the volume. Now rewrite the equality as dE = TdS — {P + E — TS — 
3/inb)(iy/V^ + 3/x(inb which implies P + E — TS — ?iiin\y = and dE/dn\^ = 3/i. 
At zero temperature, the entropy is zero (the third law of thermodynamics), 
these two equalities become 

P^-E + 3/xnb, dn^ = dE/{3i^). (6) 



In the present approach, the energy density is given in Eq. (3) from which one 
has dE — ^dpf + ^dm. Substituting this into the second equality of Eq. (6) 
and then integrating over both sides will give Eq. (5) naturally. 

Defining Ei = (i?i)nb ~ {Hi)o, the interacting part of the equivalent mass, 
mi in Eq. (2), can be re-written as mi = Ei/{3n*)/{l — {qq)rn,/ {qq)o) with 
n* = {2/3)\qq\o. Solving for the ratio {qq)nb/ {qq)o, this equation leads to 

(gg)«b ^ ^ _ El 
{qq)o 3n*mi' 



According to the expression of i^qcD, the total energy density of the quark 



4 



system can be expressed as 



E= " 



^ J \Ip' + mop'dp + Ei = Eo + Ei. (8) 





The first term is the energy density without interactions, the second term Ei 
is the contribution from interactions, and the Fermi momentum for the 
non-interacting case is connected to density by pfo = [(187r^/5')nb]^^"^ . 

On the other hand, E has already been expressed in Eq. (3). So, replacing the 
E on the left hand side of Eq. (8) with the right hand side of Eq. (3), then 
dividing by 3nb, one has 



,Pfo/ \mj \moJ 3nt 



Substituting Eq. (9) into Eq. (7) leads to 



(9) 



1% = 1 - ^ [4-^^ (-) - ^oF 
{qq)o n*m\ [pjQ \mJ \mo 



(10) 



At the same time, taking the derivative with respect to nb on both sides of 
Eq. (3) gives 



Therefore, if one knows the energy density E. the equivalent mass m and 
effective Fermi momentum pf can be obtained by solving the Eqs. (11) and 
(3), and the chiral condensate can then be calculated from Eq. (10). The great 
advantage of this scheme to calculate the in-medium chiral condensate is that 
one does not need to make any assumption on the current mass derivatives of 
model parameters. In fact, the current mass dependence can be derived like 
this. Substituting Eq. (3) into Eq. (1), performing the derivative with respect 
to m at fixed pf, then comparing with Eq. (10), one will have 

dm _ mF{pf/m) - (pfo/pf)^moF(pfo/mo) 
drrio 'mif{Pi/Tn) 

with f{x) = llxVx'^ + 1 - \n{x + \)\lx^. 
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Let's compare the present equation of state to the famous MIT bag modeL 
According to Eqs. (6), (3), and (5), the pressure of a cold quark plasma is 



487r^ m V 




- 1 (2^ -5]+ 3ch 

\ m"^ J 

- m? ch~^ ( — ) 
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1 / 



m 
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mdm, 



(13) 



where ch~^(a;) = ln(a; + ^/x"^ — 1). For comparison purpose, the independent 
state variable has been changed to the chemical potential n through Eq. (4). 
Obviously, the second term is the correspondence of the bag constant B. The 
difference is that the bag constant here is not really a constant. Instead, it is 
density-dependent . 

Because baryon masses depend on the in-medium chiral condensate, the total 
energy of the system also depends on it [8]. In principle, the hadronic matter 
bellow chiral phase transition point, where the quark condensate is nonzero, 
should be described in terms of hadronic degree of freedom. Otherwise, both 
perturbative and nonperturbative quark interactions, including quark confine- 
ment, should be accounted [9]. So let's expand the equivalent mass to Laurent 
series (mj must be divergent at pf = or nb = due to quark confinement), 
and merely take the leading terms in both directions: 

mi^ — + aipf. (14) 



In the following, it will be seen that the first term is non-perturbative, mainly 
originated from the linear confinement of quarks, and the second term is from 
the leading contribution of perturbative interactions with the coefficient Oi 
connected to the QCD coupling ag by 



It is not possible to compare the present formalism at finite density with 
lattice results. However, there are several expressions for the pressure of a cold 
quark plasma, e.g., those from the hard-thermal-loop perturbation theory [10] 
and from the weak-coupling expansion [11,12]. Although they are different 
in higher orders, their leading term is identical. Let's assume the interacting 
equivalent mass mi is, at the perturbative densities and to leading order, also 
proportional to the chemical potential fi, i.e., mi = aQfi. It is shown in the 



following that the coefficient is ao = ^J2as/n. In fact, the pressure in Eq. (13) 
can be expanded to Taylor series at mi = 0, i.e.. 




(15) 
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where Pid = g / {AST:'^)[fi^J fi^ — mo (2|U^ — bml) + Smgch ^(/^/mo)] is the pres- 
sure of a degenerate non-interacting quark plasma. The convergence of Eq. 
(16) can be mathematically proven. Now substituting mi = a^ii = ^j2a^/niJ, 
into the above perturbative expansion, and then taking the limit of mo 
due to the extreme smallness of the current mass of light quarks, the sec- 
ond term vanishes while the first and third terms give P/Pid = 1 — ^as/ir 
which is consistent with the hard-thermal-loop resumed pressure [10] and the 
weak-coupling expansion [11,12]. 

At the same time, the second term of Eq. (14) dominates at higher densities, 
i.e., mi = aipf. So inserting pf = mi/ai into Eq. (4) and solving for mi give 

mi = ai[ 
of mo 



1 + af )/i^ — ml — aimo]/ (1 + af) which approaches, at the limit 
0, to mi 



aifi/ ^1 + af. Comparing this with mi = apfi gives 
+ al = ao- Solving for oi from this equality then leads to Eq. (15). 

In order to perform the integration in Eq. (5), one should know the total 
derivative of Eq. (14), i.e.. 



dmi dmi dmi dai 
dpf dpf dai dpf 



a_i dai 
— Y + ai +pf 



Pf 



dpf 
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From Eqs. (4), (15), and the Gell-Mann-Low equation iidas/d/i = J2'^=o CnC^g'^^ 
[13], it can be shown that the derivative of ai with respect to pf is 



dai 
dpf 



Pf + mdm / dpi 



das /dai / o _i_ ™2^ _ rn ^'^ 

^^-^-^(Pf+mj m^ 



(18) 



with dm/dpf — —a-i/pj -\- oi, dm/dai — pf, da^/dai — 7rai/(l -|- af)^. To 
second order, the Gell-Mann-Low equation can be integrated out, giving 

77^ + TT^'^—F, TT-i hColn— = 0, (19) 

a, a,(l) Co Ci + Co/a, ^ A ' ^ ' 
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where as(l) is the value of ag at the QCD scale point // = A, Co = — (llA^c ~ 
2Nf)/{QTT) and Ci = -[(34A/'c- 13iVf)iVe + 3iVf/Ar,]/(247r2) arc the Gell-Mann- 
Low coefficients [14]. The scale parameter A is usually taken to be 300 MeV 
while as{l) is taken to be 1 [11]. The confinement parameter a_i should satisfy 
the constraint: the energy per baryon for the two-flavor case is no less than 930 
MeV, in order not to contradict standard nuclear physics [15]. In the present 
calculation, a_i is taken to be such a vahic that the maximum value of the 
QCD running coupling ag just does not exceed 1 at whole densities. If larger 
a_i values are used, the condensate goes to zero more rapidly. 

For a given pi, the running coupling as is obtained by solving Eq. (19) with 
the help of Eqs. (4), (14), and (15). Inserting Eq. (17) with Eq. (18) into Eq. 
(5) will give the baryon number density nt. The chiral condensate can then 
be calculated from Eq. (10). 

Now let's find the lower density behavior of Eq. (10). Because m f« mj ^ 
Pi and F{p{/m) ^ 1 at lower densities, Eqs. (11) and (9) become mj = 
dEi/{3dni,) and {pi/p{o)^mi = Ei/^Stiy,). Denoting the inter-quark interac- 

1 /3 

tion by v(f) with f cx l/n^, being the average distance of quarks [3]. The 
interacting energy per baryon Ei/{3ni,) is proportional to v(f). Because the 
confinement interaction dominates at lower densities while the confinement 
is linear according to lattice simulations [16], one has -Ei/(3nb) oc f^. Here 
the confinement exponent is denoted by Z. For linear confinement, it is equal 
to unity. With a proportion coefficient a, one can write Ei — 3crn^ . Sub- 
stituting this into the approximate equalities of Eqs. (11) and (9) just ob- 
tained gives mi = a{l - Z/3)/n^^^ and pf = [(187rV^)nb/(l - Z/S)]^/^. These 
two equalities give an accoimt for the validity of the first term in Eq. (14). 
They also show that the effective Fermi momentum has been boosted to a 
higher value by a factor, i.e., Pf/pfQ = 1/(1 — ^/3). This ratio can also be 
obtained by expanding the integrand of the integration in Eq. (5) to Taylor 
series at pf = 0, taking the leading term 2pf/(3m), making the variable sub- 
stitution mi = a^i/pf [Z=l corresponds to the first term of Eq. (14)], and 
then completing the integration. Finally substituting the ratio into Eq. (10) 
gives {qq)n^/{qq)o = 1 - Wb/ni with m = (1 - Z/2>)n*. 

Now demanding a compatibility with the famous model-independent result in 
nuclear matter, i.e., {qq)p/ {qq)o = ^ ~ p/p* with p* = M^F^/a-^ [2], one can 
get, from ni = p*, an interesting relation 

--1^ (20) 



which relates the pion-nucleon sigma term a^, the confinement exponent Z, 
and the average current mass oi u/d quarks whose value is about mo = 7.5 
MeV. In obtaining the above equation, the Gell-Mann-Oakes-Renner relation 
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~2mo{qq)Q — M^F^ has been applied. For the linear confinement of Z — 1, 
Eq. (20) gives = ^trq ~ 34 McV. The reported o"n values in literature are 
significantly different. Previously, it was reported to be a smaller one, e.g., 
in the range of 25 — 26 MeV [17]. Later, a bigger value, e.g., un ~ 56.9 ± 
6.0 MeV [18] was arrived at. After that, modest values of about 45 MeV 
was obtained [19,20]. The accepted value used to be as high as around 65 
MeV [21]. Recently, a different value of about 50 MeV was reported again 
[22]. This phenomena is due to the fact that ctn is not directly measurable 
but extrapolated according to special models. Obviously, the (Tn value here 
is bigger than previously determined, but smaller than presently accepted. 
Naturally, if one would like to reproduce the bigger values, the confinement 
would have to be also bigger. However, lattice simulations favor the linear 
confinement [16]. 
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Fig. 1. (a) Density dependence of the in- medium chiral condensate. It is generally a 
decreasing function of the density (star line), approaching to zero at about 4 fm"^. 
The open-circle Une on the left is the linear extrapolation. The up-most line with 
crosses is for a quark system without interactions. The plus-marked line gives the 
corresponding equivalent mass, in unit of the nucleon mass 939 MeV. The star line 
is the QCD coupling Og. (b) Velocity of sound in the pure perturbative calculation 
(the dotted line) and in the present approach (the dashed line). 

Numerical results are plotted in Fig. 1. In Fig. 1 (a), the up-most line is 
for the chiral condensate in a free quark gas. Because the decreasing speed is 
proportional to the current mass which is very small, the line goes down slowly. 
When interactions are considered, the condensate (full circles) decreases much 
more rapidly, approaching to zero at about 4 fm~^ which is still an extremely 
high density. The hnear behavior is shown with an open-circle hne on the left. 
The equivalent mass is given, in unit of the nucleon mass Mn = 939 MeV, 
with a plus-marked line, which shows that the ratio m/M^ is about 0.3-0.4. 
In constituent quark models [23] , the constituent quark mass is usually taken 
to be Mn/3, consistent with the present result. 
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Equations (13) and (3) provide the equation of state. If the energy per baryon 
-E/nb is plotted as a function of the density rZb, the minimum point should 
correspond exactly to the zero pressure, which is a general requirement of 
thermodynamics, as has be shown in Ref. [24]. However, this is not the case 
in the pure perturbative results, e.g., a check of the figure 7 in Ref. [11] shows 
that the density corresponding to zero pressure is smaller than the density cor- 
responding to the minimum energy per quark. This has a strong consequence 
on the velocity of sound, calculated by \dPjdE\^^'^^ which has been shown in 
Fig. 1 (b) for the present approach (the dash line) and in the weak expansion 
(the dotted hue). At higher densities, they approach asymptotically to the 
ultra-relativistic case (the full line) as expected. However, their lower density 
behavior is completely opposite. This may have some stringent consequences 
on dynamic situations, such as the possibility of strange candidates in neutron 
stars or perhaps heavy ion collisions. 

In summary, the density behavior of the in-medium chiral condensate has been 
studied with a new approach. The chiral condensate is shown to be generally a 
decreasing function of density. At lower densities, it decreases linearly, which 
means that the pion-nucleon sigma term is 9/2 times the average current 
mass of hght quarks if the confinement is hnear. With increasing densities, 
the deviation from linear extrapolation becomes significant. Considering both 
perturbative and non-perturbativre effects, the condensate approaches to zero 
at about 4 fm~^. 
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